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1  Summary 

Providing  pilots  with  a  heads  up  display  (HUD)  capable  of  presenting  a  3-Dimensional  (3-D) 
scene  to  a  pilot  during  instrument  landing  would  be  beneficial  for  spontaneous  night  landings  in 
unknown  locations.  Remote  sensing  via  RADAR  has  the  opportunity  to  provide  data  that  fulfdls 
this  capability.  The  HUD  would  be  used  to  display  information  about  barriers  and  obscurations 
in  the  flight  path  and  intended  landing  area. 

This  report  covers  the  transformation  and  sensor  requirements  necessary  for  RADAR  imagery  to 
appear  in  the  same  geometry  as  an  equivalent  optical  image;  what  the  pilot  sees.  Spatial 
resolution,  slant-range  to  ground-range  transformations,  and  maximum  detectable  range  are 
covered.  Appendices  are  provided  which  contain  more  information  on  interferometric  synthetic 
aperture  RADAR  (IFSAR)  and  an  approach  for  removing  terrain  artifacts  from  the  data 
transformations. 


2  Introduction 

The  overall  goal  of  the  “RADAR  Imaging  Transformation  for  Heads  Up  Display”  study  is  to 
explore  ways  in  which  a  pilot  could  perform  an  instrument  only  landing  using  RADAR  images 
as  the  main  source  of  information.  RADAR  can  be  used  in  most  types  of  weather  and  during 
night  landings,  making  it  a  logical  choice  for  spontaneous  landing  of  aircraft.  The  RADAR  data 
must  be  displayed  in  a  manner  that  is  consistent  with  how  a  pilot  looks  out  of  a  cockpit.  The 
geometry  of  the  display  must  be  the  same  as  the  geometry  of  the  real  world  through  the  human 
visual  system. 

RADAR  data  collections  take  place  in  a  range-angle  geometry,  as  RADAR  is  typically  used  as 
measure  of  distance.  Optical  images  are  collected  using  an  angle-angle  geometry,  and  therefore  a 
transformation  is  needed  to  convert  RADAR  imagery  into  its  equivalent  optical  image.  This 
report  serves  to  defining  sensor  requirements,  relate  electro-optic  sensor  parameters  to  the 
equivalent  RADAR  parameters,  and  formulate  the  data  transformations  for  appropriate  display  of 
RADAR  imagery. 


3  Assumptions 

3. 1  Flat  Earth  Assumption  in  RADAR 

The  flat  Earth  assumption  used  in  RADAR  is  made  because  normally  the  radius  of  Earth  is  much 
larger  than  the  distance  to  the  target.  The  assumption  states  that  the  affect  of  the  spherical  shape 
of  Earth  is  insignificant  when  one  is  close  to  the  surface  of  Earth.  The  parameters  below 
describe  the  landing  scenario  used  throughout  this  study,  with  Figure  1  shown  to  help  guide  the 
derivation  of  the  measurement  error.  The  difference  in  straight-line  distance  to  the  runway, 

Rmax  meas  or  Rr,  and  arc  length  along  the  curvature  of  the  Earth,  Rarc,  determines  the  possible  error 
in  measurement. 
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Figure  1:  The  geometry  used  for  deriving  the  error  inherent  in  using  the  flat  earth  assumption 


Based  on  the  above  geometry,  the  angle,  <9,  can  be  calculated  as  shown  in  equation  (1).  Once 
6>is  known,  arc  length  can  be  calculated  using  equation  (2). 
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Given  the  radius  of  Earth,  RE,  as  6371  km  and  the  range  to  the  runway,  Rr,  as  1.635  km,  the 
measurement  error  can  be  calculated  to  be  4.48  jam,  as  shown  below. 

0  =  0.0147° 

0  =  2.56  xlO^4  rad 


Rare  —  ^  RE 

Rare  =  1 -63500000448km 
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Rarc  =1.63500000  448km 

Rr  =  1.635km 

AR  =  Rarc  -  Rr  =  4.48  pm 


As  can  be  seen  from  the  above  calculation,  for  the  parameters  given,  the  flat  Earth  assumption  is 
valid.  The  difference  in  the  range  to  the  runway  between  straight  line  distance  and  arc  length  is 
4.5um.  Note  that  this  discussion  pertains  to  the  curvature  of  Earth,  and  not  to  terrain  models. 
Terrain  models  are  discussed  in  Appendix  C. 


4  RADAR  and  EO  Sensor  Similarities  and  Differences 

The  following  discussion  covers  the  differences  between  electro-optics  (EO)  sensors  and 
RADAR  sensors,  specifically  in  the  areas  of  data  collection  style  and  spatial  resolution.  In  order 
to  better  transform  the  RADAR  image  to  its  optical  equivalent,  height  information  needs  to  be 
taken  into  account.  Interferometric  synthetic  aperture  RADAR  (IFSAR),  can  be  used  to  estimate 
the  heights  of  objects  within  a  scene  and  is  discussed.  The  maximum  detectable  range  of  a 
RADAR  system  is  also  explored. 


4. 1  RADAR  Mode  Definitions 

The  transformation  from  a  RADAR  geometry  to  an  EO  geometry  can  be  described  in  terms  of  a 
B-scope  to  C-scope  conversion.  The  B-scope,  as  shown  in  Figure  2,  displays  an  image  in  range 
and  azimuth,  the  native  form  of  RADAR  imagery.  A  C-scope  image  has  undergone  a 
perspective  transformation  and  appears  as  an  EO  image  would  appear  to  a  pilot.  The  difficulty 
in  this  transformation  comes  in  correctly  relaying  height  information  in  the  C-scope  image. 


B  scope 


Azimuth 

Figure  2:  An  example  of  B-scope  imagery  and  C-scope  imagery 
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4.2  Collection  Differences 

The  collection  geometry  of  RADAR  systems  differs  from  the  geometry  of  electro-optical 
systems.  Electro-optics  data  is  collected  such  that  pixels  are  arranged  on  a  Cartesian  grid.  Pixel 
location  within  the  grid  is  determined  by  noting  two  angles;  azimuth  angle  ai  and  depression 
angle  012,  as  shown  in  Figure  3,  in  order  to  locate  the  point  along  line  L/.  Due  to  this  setup, 
optical  systems  can  also  be  referred  to  as  angle-angle  systems.  The  disadvantage  of  an  angle- 
angle  system  is  that  an  object  must  be  observed  from  two  positions  in  order  to  image  in  three 
dimensions. 

Pixel  size  and  quantization  determine  the  image  quality  and  are  dependant  upon  the  imaging 
hardware. 


RADAR  systems  create  images  by  calculating  the  range  to  an  object  and  noting  the  angle  in 
azimuth  and  depression  of  the  receiver  relative  to  a  fixed  reference  point.  This  range-angle 
geometry  allows  for  three  dimensional  imaging  from  a  single  observation  point.  Range  bins  exist 
as  concentric  circles  about  the  receiver,  Due  to  this  geometry,  interpolation  is  necessary  in  order 
to  view  the  RADAR  pixels  correctly  on  a  Cartesian  grid.  Figure  4  demonstrates  the  collection 
geometry  of  a  RADAR  system,  in  this  case  a  system  used  for  synthetic  aperture  RADAR  (SAR). 
Point  Pi  is  identified  by  its  location  in  the  circular  range  bin  C/  by  the  range  R,  and  a(i,  the  angle 
between  the  sensor  path  and  Pi. 
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A  comparison  of  the  angle-angle  geometry  of  EO  and  the  range-angle  geometry  of  RADAR  can 
be  seen  in  Figure  5.  An  EO  image  provides  intensity  information  about  a  single  plane  within  the 
scene,  as  each  point  in  the  plane  can  be  referenced  by  its  depression  angle  and  azimuth  angle. 
Information  about  the  distance  of  a  point  from  the  reference  point  cannot  be  directly  obtained. 
RADAR  systems  inherently  detect  range,  therefore  a  point’s  location  along  the  dotted  lines  in 
Figure  5  is  provided  along  with  the  point’s  coordinates  in  azimuth  and  depression  angle. 
Therefore,  the  advantage  of  the  range-angle  geometry  is  the  ability  to  obtain  three-dimensional 
information  (depression  angle,  azimuth  angle,  and  range)  from  a  single  reference  point. 


Figure  5:  RADAR  imagery  provides  an  extra  dimension  of  information  over  EO  imagery 
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4.3  Spatial  Resolution 

Spatial  resolution  determines  how  well  one  can  distinguish  features  in  an  image.  In  both  optics 
and  scanning  RADAR,  spatial  resolution  is  dependant  on  frequency  and  the  dimensions  of  the 
equipment  used  for  signal  detection. 

The  spatial  resolution  for  optics  shown  in  equation  (3)  relies  upon  Rayleigh  resolution  criteria 
with  the  diffraction  of  light  as  the  limiting  factor.  This  equation  assumes  a  spherical  lens  system; 
therefore  the  resolution  in  both  of  the  cross-range  dimensions  is  the  same. 


p  =  1.22 — 

D 


A  =  wavelength  (  3  ) 

d;  =  distance  from  lens  to  image 
D  =  diameter  of  aperture 


Due  to  the  nature  of  RADAR  collections,  spatial  resolution  is  specified  in  both  range  and  cross 
range.  Cross  range  resolution  in  RADAR,  as  shown  in  equation  (4),  is  similar  to  the  optical 
resolution  of  equation  (3).  Both  equations  are  dependant  on  wavelength,  distance  from  sensor  to 
target,  and  the  main  dimension  of  the  detecting  hardware. 

p.,=m=^R 

P  =  beamwidth  (  4  ) 

A  =  wavelength 

R  =  range 

D  =  antenna  dimension 


Range  resolution  in  scanning  RADAR  is  a  function  of  the  pulse  length  of  the  RADAR  sensor,  as 
shown  in  equation  5.  Pulse  length  and  bandwidth  have  an  inverse  relationship. 


c 

2B 


c  =  speed  of  light  (  5  ) 

t  =  pulse  length 
B  =  bandwidth 


Another  form  of  RADAR  imaging,  SAR,  can  provide  increased  azimuth  resolution  by  combining 
several  collections  from  a  moving  platform  into  a  single  beam,  effectively  increasing  the  antenna 
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dimension  and  creating  a  finer  azimuth  resolution.  Figure  6  shows  a  comparison  of  the  resolution 
cell  size  for  scanning  RADAR  and  SAR.  The  synthetic  beamwidth  in  the  SAR  case  is  a  result  of 
the  processing  of  several  data  collects. 


Scan  RADAR 


SAR 


Compressed 
pulse  length 

Resolution 

cell 


Synthetic 

beamwidth 


Compressed 
.'l  impulse  length 

Resolution 
cell 


Figure  6:  SAR  processing  allows  for  better  resolution  than  scanning  RADAR 


Equation  (6)  and  equation  (7)  show  the  updated  spatial  resolution  equations  when  a  SAR 
processing  is  used.  The  mainlobe-broadening  factors,  Ka  and  Kr,  are  constants  associated  with 
the  antenna  hardware. 


Paz 


2A  9 


Ac  =  center  frequency 
Ka  =  mainlobe  -  broadening  factor 
A0  =  coherent  integration 
angular  interval 


(6) 


Pr  = 


c_K^ 

2  B 


c  =  speed  of  light  (  7  ) 

Kr  =  mainlobe  -  broadening  factor 
B  =  pulse  bandwidth 
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4. 4  Spatial  Frequency  Relationship 

Another  way  of  interpreting  the  difference  between  EO  and  RADAR  is  in  the  spatial  frequency 
domain  of  the  scene.  In  the  spatial  frequency  domain,  an  optical  image  occupies  a  plane  parallel 
to  the  fx-fy  plane  shown  in  Figure  7.  Due  to  the  differences  in  collection  geometries  stated 
earlier,  SAR  data  collections  fill  in  the  three  dimensional  space  by  collecting  data  around  a  cone 
as  shown  in  Figure  7.  A  two-dimensional  cut  through  this  cone  shows  that  a  RADAR  sensor  will 
fill  in  only  one  line  in  the  same  plane  as  an  optical  sensor. 


f. 


y 


Figure  7:  Spatial  frequency  relationship  between  an  optical  sensor  and  SAR 


4.5  Interferometric  Synthetic  Aperture  RADAR  (IFSAR) 

Interferometric  synthetic  aperture  RADAR  (IFSAR)  requires  the  use  of  two  receivers  offset  from 
one  another  by  a  known  amount.  This  collection  geometry  allows  for  the  discernment  of  height 
information  based  on  the  phase  differences  between  the  received  signals  on  the  two  antennas 
from  a  single  pulse  aimed  at  the  same  target.  A  derivation  of  the  height  equation,  shown  in 
equation  (8),  can  be  found  in  Appendix  A.  In  equation  (8),  z  is  the  height  of  an  object,  hant  is  the 
height  of  the  antennas,  and  R  is  the  slant  range  between  an  antenna  and  the  target.  The  imaging 
geometry  shown  in  Figure  8  is  representative  of  the  collection  geometry  for  the  example  data  set 
used  in  algorithm  testing  for  the  overall  project. 


z  =  ha,u  -RcosO  (m) 
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(8) 


hant  =  antenna  height  =  100  m 

Transmit  from  one  antenna  B  =  interferometer  baseline  =  0.024  m 

and  receive  on  both  antennas 

p  =  range  to  target  =  [0... 3615.2]  m 


two  antennas 

Figure  8:  Data  Collection  Geometry  used  in  this  study  for  algorithm  development 


4. 6  Slant  Range  to  Ground  Range  Transformation 

Slant  range  to  ground  range  distortion  in  RADAR  imagery  occurs  because  the  RADAR  is 
measuring  distance  along  the  slant-range  as  opposed  to  the  true  horizontal  ground  distance  to  a 
target.  As  shown  in  the  right  side  of  Figure  9,  this  causes  targets  in  the  near  range  to  appear 
compressed  relative  to  the  far  range. 

A  projection  operator  such  as  the  one  in  equation  (9)  can  be  used  to  convert  slant  plane  imagery 
into  ground  plane  imagery. 


Figure  9:  Slant  range  to  ground  range  conversion  geometry  with  associated  imaging  artifacts 
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4. 7  Maximum  Detectable  Sensor  Range 

The  hardware  of  a  RADAR  system  determines  the  maximum  detectable  range  for  a  given  target 
signal-to-noise  ratio  (SNR).  This  maximum  detectable  range  is  analogous  to  the  horizon  of  an 
optical  image,  which  is  limited  by  the  curvature  of  Earth.  Because  a  RADAR  sensor  is  its  own 
source,  the  RADAR  horizon  can  be  defined  as  the  minimum  of  the  detectable  range  due  to 
geometry  and  detectable  range  due  to  the  output  power  of  the  RADAR  transmitter  for  a  given 
SNR. 

Appendix  B  contains  a  short  primer  on  the  basic  RADAR  equations  necessary  to  construct  an 
estimate  for  the  maximum  detectable  sensor  range  based  on  parameters  used  in  this  study.  The 
derivation  for  this  equation  also  necessitates  a  description  of  the  RADAR  cross  section  of  a 
planar  target,  which  is  contained  in  section  4.7.2. 


4.7.1  Maximum  Detectable  Range  Due  to  Collection  Geometry 

A  good  estimate  for  the  geometric  RADAR  horizon  is  provided  in  the  Naval  Electronic  Warfare 
and  RADAR  Systems  Engineering  Handbook  and  is  shown  in  equation  (10).  The  value  %  is  an 
empirical  value  used  to  estimate  the  wave  distortion  of  the  atmosphere,  RE  is  the  radius  of  Earth, 
and  hant  is  the  height  of  the  antenna. 

R„„  =  V2-X'Vh-  0»)  (10) 


4.7.2  RADAR  Cross  Section  for  Planar  Target 

The  RADAR  cross  section  (RCS)  of  an  object  is  a  measure  of  how  an  object  reflects  an  incident 
electromagnetic  wave.  When  trying  to  calculate  the  distance  to  the  RADAR  horizon,  there  is  not 
one  specific  object  reflecting  energy,  as  the  Earth’s  surface  as  a  whole  must  be  taken  into 
account.  This  can  be  done  if  the  Earth’s  surface  is  treated  as  distributed  clutter  over  a  planar 
target.  The  RCS,  oh  of  the  Earth’s  surface  is  defined  by  the  back-scatter  coefficient  of 
distributed  ground  clutter  and  the  resolution  cell  of  the  RADAR  as  shown  in  equation  (11). 

Figure  10  shows  the  geometry  for  this  scenario. 

<*t  =/VA,z-°bs(m2)  (11) 

pr  =  range  resolution  (m) 
paz  =  azimuth  resolution  (m) 
abs  =  back  -  scatter  coeff.  (unitless) 

Px  =  beamwidth  (rad) 

R  =  range  (m) 

B  =  bandwidth  (Hz) 


Figure  10:  Geometry  for  calculating  the  RCS  of  a  planar  target 
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The  backscatter  coefficient  of  distributed  clutter  depends  on  the  angle  of  incidence  and  has  been 
previously  measured  by  Gatesman,  et  al,  as  shown  in  Figure  11. 


Figure  11:  Back  scattering  coefficient  for  distributed  ground  clutter 


4.7.3  Maximum  Detectable  Range  Due  to  Transmitter  Power 

RADAR  range  depends  on  transmitter  power  because  a  RADAR  is  its  own  source.  The 
maximum  range  for  a  given  transmit  power  can  be  found  by  solving  the  standard  RADAR 
equation,  equation  (12)  for  range. 


PtGtGrX2Gpqt 

(4ti)3R2R2 


(12) 


Assuming  a  monostatic  collection  geometry  and  replacing  at  with  equation  (11),  equation  (12) 
becomes: 


P  = 


rg2x2g„  V 


(4tt)3R4 


\ 


v55'|iR'‘\ 


(13) 


Using  the  antenna  gain  definition  from  Appendix  B  and  solving  equation  (13)  for  R  provides  the 
equation  for  the  maximum  detectable  RADAR  range  in  terms  of  transmit  power  in  equation  (14). 


R 


max 


Ptdxd^2fcqbsGp 

8;t  B  •  SNR  •  Pn 


(14) 
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4.7.4  Calculation  of  RADAR  Horizon  using  given  values 

The  set  of  parameters  in  Figure  12  was  used  in  this  study  for  algorithm  development. 


hant=100m 
dx  =dy  =0.12m 

=  0.85  (typical  estimate) 
fc  =  94  GHz 

abs  =  -8  dB  =  0. 1 6  (unitless) 

Gp  =  16384  =  42  dB  (unitless) 

B  =  300  MHz 

SNR  =  [20, 35, 45  dB]=  [100,3162,31623] 

Pn  =-llldB  =  7.8xlO’12W 

Figure  12  :  Given  parameters  for  the  RADAR  collection  used  in  this  study 

Using  these  parameters  in  equation  (10)  gives  a  maximum  range  due  to  geometry  of  41,260 
meters.  Solving  equation  (14)  for  the  maximum  range  due  to  transmit  power  gives  a  distance  of 
3722  meters,  assuming  a  transmit  power  of  one  Watt  and  an  SNR  of  20  dB.  Therefore,  for  this 
scenario,  the  RADAR  transmit  power  is  the  dominant  contributor  to  the  overall  maximum 
detectable  range.  Figure  13  provides  the  maximum  detectable  sensor  range  for  this  experiment 
for  three  values  of  SNR  and  a  range  of  transmit  powers. 


Maximum  Detectable  Range 


Figure  13  :  Maximum  detectable  range  for  the  scenario  in  this  study 
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5  Conclusions 

The  purpose  of  this  report  was  to  examine  some  of  the  required  RADAR  sensor  parameters  and 
transformations  in  order  to  convert  a  RADAR  image  into  its  equivalent  optical  image  for  the 
purpose  of  landing  an  aircraft  via  a  heads  up  display.  Equations  for  spatial  resolution,  slant 
range  to  ground  range  transformation,  and  maximum  sensing  distance  were  presented. 

Concern  over  adjusting  for  the  underlying  terrain  of  an  image  could  be  addressed  as  shown  in 
Appendix  C.  Future  work  on  this  project  could  involve  furthering  the  ideas  of  Appendix  C  to 
create  an  automated  way  of  adjusting  for  the  terrain  in  an  image. 
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Appendix  A  -  Derivation  of  Height  Equation  Using  IFSAR 


Each  pixel  in  an  image  contains  two  data  points,  one  from  each  of  the  receive  antennas.  The  two 
received  signals  are  shown  in  equation  (Al)  and  equation  (A2),  with  the  major  difference  in  the 
signals  being  the  extra  phase  term  in  I2.  The  geometry  for  this  derivation  is  the  same  as  in  Figure 
8. 


/[  =  axe 


2  71 

J'\  <Pri+j^2r 


(Al) 


I2  =  a2e 


j\  Vri* — — 


(A2) 


If  the  interferometer  spacing  B  is  not  too  large,  the  amplitude  and  phase  of  the  complex 
reflectivity  is  nearly  the  same  for  each  antenna,  thus  the  simplifications  a  =  at  =  a2  and 
<pR i  =  (pR2  can  be  made. 

If  we  form  an  interferogram  from  the  images,  we  have: 

A, 2  =  I A 

=  a2e[  A  J  (A3) 

/l2  =  a2ei<p 

,  2nAR 

where  (p  = - . 

2 


The  law  of  cosines,  equation  (A4),  allows  us  to  obtain  AR  by  expanding  and  dropping  terms 

2  ( B2'\ 

0(AR  ),  O  —  .  The  resulting  equations  for  AR,  and  ultimately  (p,  are  shown  in  equation  (A5) 


yRj 


and  equation  (A6)  respectively. 


(R  -  A Rj  =R2+B2-  2  RB  cos 

AR  =  B  sin(#  -  a) 
2xB  . 

(p  = - sint^-o:) 

A 


n 


—  0 + oc 


( A4 ) 
(A5) 
(  A6 ) 


Solving  for  6  shows  that  if  we  measure  the  phase  at  each  image  point  we  can  obtain  the  angular 
location  of  the  image  points.  This  relationship  is  explicitly  shown  in  equation  (A7). 
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6  =  sin 


-i 


A  (p 

2  7rB 


+  a 


( A7  ) 


We  can  obtain  the  y  and  z  coordinates  of  each  image  point  using  equation  (A8)  and  equation 
(A9). 

y  =  R  sin  6  ( A8  ) 

z  =  H  —  R  cos  0  (A9) 

In  summary,  if  a ,  B,  and,  H  are  known,  then  height  of  an  object,  can  be  calculated  from  the 

interferogram,  as  height  is  dependant  upon  6,  which  is  derived  from  cp,  which  is  determined  by 
the  interferogram. 
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Appendix  B  -  Basic  RADAR  Equations 


Noise  Power 

P  =P  =KTFB 

n  noise 

K  =  1.38x1  O'23  W/(Hz-K) 

T  =  system  temperature  (K) 

F  =  loss  factor  (unitless,  ~  8  dB) 
B  =  bandwidth  (Hz) 

Antenna  Beamwidth 


(3  =  beamwidth  (rad) 

X  =  wavelength  (m) 
d  =  antenna  dimension  (m) 

Antenna  Gain 

l2 

dx  =  antenna  dimension  (m) 
dy  =  antenna  dimension  (m) 
t;  =  antenna  efficiency  (unitless) 

Signal-to-Noise  Ratio  (SNR) 

SNR  =  — 

Pn 

SNR  =  signal  to  noise  ratio 
Pr  =  received  power(W) 

Pn  =  noise  floor  (W) 


(Bl) 


(  B2 ) 


(  B3  ) 


(  B4 ) 
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RADAR  Equation 


Figure  B-l:  RADAR  Equation 


P,G,G,X;Gpg, 

'  M’rX 

Pr  =  receiver  power  (W) 

Pt  =  transmitter  power  (W) 

Gt  =  antenna  gain  (unitless) 

Gr  =  antenna  gain  (unitless) 

X  =  wavelength  (m) 

Gp  =  range  compression  processing  gain  (unitless) 

ot  =  target  RADAR  cross  section  (RCS)(m2) 

Rt  =  range  from  antenna  to  target  (m) 

Rr  =  range  from  target  to  antenna  (m) 


(B5) 
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Appendix  C  -  Algorithm  for  Testing  Terrain  Distortion 


A  possible  extension  of  this  program  is  to  include  terrain  modeling  into  the  RADAR  imagery  in 
order  to  better  mimic  the  optical  imagery.  As  shown  in  the  SAR  geometry,  a  transformation 
from  slant  range  to  ground  range  is  performed  in  order  to  better  align  an  image  with  its  optical 
counterpart.  However,  this  transformation  assumes  that  the  projection  is  taking  place  onto  a  flat 
ground,  which  is  not  true  in  most  cases.  A  potential  way  to  overcome  this  problem  is  shown  in 
Figure  Cl. 


Figure  C-l:  Overview  of  algorithm  for  incorporating  local  terrain  modeling 


A  Digital  Elevation  Model  (DEM),  such  as  the  one  in  Figure  C2,  can  be  used  to  separate  hilly 
areas  of  a  scene  from  relatively  flat  areas  of  a  scene.  The  threshold  for  this  separation  has  yet  to 
be  determined. 
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Figure  C-2:  An  example  of  a  Digital  Elevation  Model  (DEM) 


The  hilly  areas  of  the  scene  can  then  be  transformed  into  a  series  of  wedges,  one  of  which  is 
shown  in  Figure  C3.  The  wedges  can  then  be  rotated  in  three  dimensions  until  a  suitable 
orientation  has  been  found  for  each  wedge. 


^flat  ^hill 

dhin’:  hill  length 

dflat’:  distance  between  tower  and  start  of  hill 
h’:  hill  height 

a’:  cross-range  location  of  hill  center 
c(>:  hill  orientation  with  respect  to  range 
0:  grade  of  the  hill 

Figure  C-3:  One  example  of  a  terrain  model  correction  wedge 
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Three  dimensional  (3-D)  rotation  algorithms  are  commonplace  in  computer  graphics  rendering 
engines.  An  example  of  the  transformation  matrices  necessary  for  computing  a  3-D  rotation  are 
shown  in  Figure  C4.  Furthermore,  if  the  object  is  not  located  at  the  origin,  the  object  must  first 
be  translated  to  the  point  of  origin,  then  rotated,  and  finally  translated  back  to  the  objects  original 
location. 

cosy  siny  0  0 

-siny  cosy  0  0 

0  0  10 

0  0  0  1 

10  0  0 

0  cos  a  sin  a  0 

0  -sin  a  cos  a  0 

0  0  0  1 

cos  p  0  -  sin  p  ( 

0  1  0  ( 

sin  p  0  cos  p  ( 

0  0  0 

Figure  C-4:  Transformation  matrices  for  performing  a  3-D  rotation 

The  accuracy  of  the  terrain  model  would  be  improved  by  minimizing  the  mean  squared  error 
between  a  feature,  e.g.,  a  road,  in  the  optical  image  and  an  object  in  the  RADAR  image.  For 
example,  after  the  road  is  first  extracted  from  the  scene  in  both  the  transformed  RADAR  image 
and  the  optical  image,  the  accuracy  of  the  terrain  model  can  be  discerned  via  the  squared  error 

between  the  angle  of  the  road  in  the  RADAR  ( a )  and  optical  images  (a). 

min 

a’,dflat,dhill,(p,(9,h' 

Finally,  iterations  over  various  terrain  model  sizes  and  orientations  can  be  done  until  a 
sufficiently  small  error  is  found. 
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